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Abstract 

Chiral tunneling through a harmonically driven potential barrier in graphene monolayer is con- 
sidered in this work. Since the quasiparticles in this system are chiral in nature, tunneling is highly 
anisotropic, we determine the transmission probabilities for the central and sidebands as the in- 
cident angle of the electron beam is changed . Furthermore, we investigate how the transmission 
probabilities change as the width, amplitude and frequency of the oscillating barrier is changed. 
An interesting result of our study is that perfect transmission for normal incidence that has been 
reported for a static barrier persists for the oscillating barrier, manifestation of Klein tunneling in 
a time harmonic potential. 
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I. INTRODUCTION 



Advancement in technology has led to active investigation of electron transport in semi- 
conductor nanostructures in time-dependent fields. The additional degree of freedom pro- 
vided by the time dependence has led to the appearance of new phenomena in electron 
transport, for a review see [l| and references therein. Engineering of the confinement po- 
tential and band structure has allowed the possibility of studying photon assisted tunneling 
(PAT), where inelastic tunneling events occur in the presence of an ac field, in various driven 
systems. This topic is not only of academic interest but also has device applications. Early 
studies of PAT include the work of Dayem and Martin who provided evidence of absorption 
and emission of photons in tunneling transport in experiments on superconducting films in 
the presence of microwave fields [2]. Soon after this, Tien and Gordon theoretically justified 
this observation [3i]. They assumed a time harmonic potential difference produced between 
the two films by a microwave field and solved the time-dependent Schrodinger equation for 
the system. Their photon assisted transport model accounted for transmission in the side 
bands in the presence of microwave radiation. The basic physical idea behind photon assisted 
tunneling is that an oscillating potential can lead to in-elastic tunneling where the electrons 
exchange energy quanta (photons) with the oscillating field. In such systems, a harmonically 
driven in time potential results in exchange of energy with electrons in the units of modula- 
tion quanta fiw, uo being the modulation frequency. Therefore, electrons at energy E can be 
transferred to the sidebands at energies E±nhw {n = 0, ±1, ±2, ..) while traversing a region 
of space subjected to such a time-harmonic potential. The prototypical tunneling structure, 
which is an essential element of nanostructures where electron tunneling is investigated, is 
a single barrier. A common model in these studies is a time-modulated potential that has a 
finite spatial profile. Standard electron transport through various types of time-oscillating 
potential regions has been studied previously. More pertinent to the work undertaken here is 
that of Buttiker and Landauer. They investigated the traversal time of particles interacting 
with a barrier with time-oscillating height 0, 0]. Furthermore, M. Wagner wrote a series 
of papers on photon assisted transport through quantum wells and barriers with oscillating 
potentials 6|]. Exchange of photons between the oscillating potential and electrons transfers 
electrons to the sidebands with a finite probability. Wagner determined these transmission 
probabilities using transfer matrix methods and discussed them as a function of the dimen- 
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sionless parameter a which is the ratio of the amphtude of the time oscillating potential 
to its modulation energy. There are other contributions to this field that are relevant to 
our work and these have been put together in|3]. Recently, single layer carbon crystals 
(graphene monolayer) were fabricated which has generated considerable interest in finding a 
material that can replace silicon in micro-electronic devices. The idea of carbon based nano- 
electronics has been around since the discovery of carbon nanotubes. The recent fabrication 
of graphene monolayer has provided another avenue for carbon based electronics. Devices 
based on photon-assisted electron tunneling require the consideration of electron transport 
in time-harmonic potentials. For graphene based PAT devices it is essential to consider 
transport of charge carriers in graphene through time-harmonic potentials. To this end, we 
undertake the study presented here realizing that quasiparticles in graphene systems are 
quite different from the standard electrons that we encounter in conventional semiconductor 
based heterostructures. At low energies, quasiparticles (electrons and holes ) in graphene are 
described by the relativistic Dirac-like equation and possess charge conjugation symmetry 
as a single equation describes both particles (electrons) and antiparticles (holes). This is 
due to the crystal structure of graphene which is a layer of carbon atoms tightly packed in 
honeycomb lattice. It can be thought of as the superposition of two equivalent triangular 
sublattices conventionally called sublattice A and B. Quantum mechanical hopping between 
these sublattices results in the formation of two cosine-like energy bands. Intersection of 
these bands near the edges of Brillouin zone (Dirac points) leads to the conical energy spec- 
trum E = ±hvFk (with the effective Fermi speed ( vp = lO^m/s). Above zero energy, the 
charge carriers in these systems are electrons which are usually termed Dirac electrons. The 
2D Dirac-like spectrum was confirmed recently by cyclotron resonance measurements and 
also by angle resolved photoelectron spectroscopy (ARPES) measurement sjs]. Recent the- 
oretical work on graphene multilayers has also shown the existence of Dirac electrons with 
a linear energy spectrum in monolayer graphene jsl. The Dirac equation implies that the 
quasiparticles in graphene are chiral, tunneling through potential barriers in these systems is 
significantly different from systems where tunneling of standard electron occurs such as the 
two-dimensional electron gas (2DEG) systems realized in semiconductor heterostructures. 
Chiral nature of particles in graphene results in quantum tunneling being highly anisotropic 
where relativistic effects such as perfect transmission through high and wide barriers can 
occur (Klein tunneling Q)Q. This occurs due to the conservation of chirality in inter- 
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action with the barrier, electrons in graphene can propagate to hole states through a high 
barrier without any damping. The study of this effect is relevant to the development of 
future graphene based devices. From a basic research point of view, graphene based sys- 
tems, due to their lower 'light speed', can be quite useful for studying relativistic effects. 
Moreover, the role of chirality can be highlighted in electron transport in graphene. In 
graphene-based systems, electronic transport through barrier structures has been recently 
investigated 11, I2I, Q, 
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i ll9|. In this work, we consider the transport 
of Dirac electrons in monolayer graphene through a harmonically oscillating in time square 
potential barrier. When standard electrons pass through a region which is subjected to 
time harmonic potentials, electronic transitions from central band to sidebands occur. Here, 
when transmission of Dirac electrons is considered, we also find transitions from the cen- 
tral to sidebands at energies E ± nhu {n = 0,±1,±2,..) and determine the transmission 
probabilities for the sidebands. Moreover, we investigate how the transmission probabilities 
change as various perimeters involved in the problem are varied with emphasis on the chiral 
nature of tunneling. 



II. FORMULATION 



We consider monolayer graphene sheet in the xy-plane. The square potential barrier 
is taken to be in the x-direction while particles are free in the y-direction. Width of the 
barrier is a, height of the barrier is oscillating sinusoidally around V with amplitude Vi and 
frequency uj. Electrons with energy E are incident from one side of the barrier in monolayer 
making an angle 0o with the x-axis and leave the barrier with energy E±nhuj (n = 0, ±1, ..) 
making angles (pn after transmission and vr — (/)„ after reflection. 

Let us consider the Hamiltonian H describing the system 

H =Ho +Hi (1) 

where Hq is the Hamiltonian for the static case where the barrier height is not changing 
with time and Hi describes the harmonic time dependence of barrier height, given by 

Ho = -ihvF(T.V + V (2) 
Hi = ViCos{ujt) (3) 
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V, V\ are the static square potential barrier and the amphtude of the oscillating potential, 
respectively. Both V and V\ are constants for < x < a with a positive and are zero 
elsewhere, a = {crj;,ay) are the Pauli matrices, vp is the Fermi velocity. 

Solutions of the Dirac equation in the absence of the oscillating potential, if q ^ = E'^, 
are given in [llj] and can be used for constructing solutions to the time- dependent problem. 
For the tunneling problem, we consider the incoming electrons to be in plane wave states 
\E'i (x, y, t) at energy E 

(x, y, t) = e'^yy ( \ ) e^'^i'^e-*^*/^ (4) 



.So e' 

where ki and ky are the x— and y— component of the electron wavevector, respectively. 
So = sgn{E) and 0o is the angle that incident electrons make with the x-axis. 

Reflected and transmitted waves have components at all energies E ± Ihu {I = 0, ±1, ..) 
since the oscillating potential barrier can give and take energy away from electrons in units 
of hhj. This change in energy causes only the x-component of momentum to change. Hence, 
wavefunctions (x, y, t) for reflected and (x, ?/, t) for transmitted electrons, respectively 
are 

(x, y, t) = e'^^y ^n[_ „ .^^ ) e-*'=i^e-'(^+''^)*/'^ (5) 

l=—oo ^ ' ' 



and 



where 



(x, y, t) = e'^^y Y^ki e*'=i"e-'(^+''^)*/'^ (6) 



(t)i = tan'^ {ky/k[) 
si = sgn{E + Itkj). 

In the barrier region, where Hi is nonzero, the eigenfunctions '^^ (x, t) of H can be ex- 
pressed in terms of the eigenfunctions \l/o(x, y) of Hq as[3.] 

{x,y,t) = Mx,y) Yl f ^ j e-^--*-^^*/^ 

n=— oo ^ ^ 

where Jn (^) is the nth order Bessel function. A linear combination of wavefunctions at 
energies E + Ihoj {I = 0, ±1, ..) has to be taken. Hence 



l=oo 



l=—oo 



Ci 



-S; e- 



le '''2^ 



(7) 



-i{n+X)u)t—iEt I h 



where 



f 



tan ^{ky/k2) 



s'l = Sgn{E + inw -V). 

The wavefunctions given in equations(jl}17]) have to be continues at the boundary. Apply- 
ing this condition at x = and x = a , i.e. (0, y, t) + (0, y, t) = (0, y, t) and 
{a,y,t) = {a,y,t) and reahzing that {e*'^'^*} are orthogonal, we obtain the following 
set of simultaneous equations: 

l=oo 



A pi<t>n _ ^ p-i4>n 

here An 



i=— oo 
/=oo 

i=— oo 



Si Jn—l 



(8) 
(9) 



Z=oo 



l=—oo 
l=oo 



he 



ikia 



+ Cie 



-ikia 



l=—oo 



huj 



(10) 



The above set has infinite number of coupled equations and contains infinite number of 
unknowns( n, I goes from — oo to oo ). This linear system of equations cannot be analytically 
solved. Nevertheless, the infinite series in these coupled equations can be truncated and a 
finite number of terms starting from — upto , where A^ > ^ , retained if we note that 
the coupling strength is determined by the quantity ^ through Bessel functions Jn {^) 
and Jn (■^)) they become negligible for order n higher than Vi jhw. Equations (PfTTi) are 
numerically solved for t„. The transmission probability for the nth sideband, T„, for which 
/c" is real and corresponds to propagating waves, is obtained from: 



^ COS0„ , ,2 

Tn = ^ tn f 12 

COS 00 

whereas imaginary A;" corresponds to evanescent waves that carry no particle current with 
the resuh T„ = 0. ki can be real or imaginary depending on the particular values of 
the following parameters: incident energy E, oscillation frequency u, incident angle 0o- 
The numerical results obtained are discussed in the next section. Furthermore, analytical 
results are obtained if we consider small values of a = and include only the first two 
sidebands at energies E ± huj alongwith the central band at energy E. Moreover, we have 
to invoke the conditions hw < E such that sgn{E ± hw) = +1 and hw < \E — V\ such that 
sgn{E — V± hu) = —1 for E < V. Hence, we are able to truncate the sums in equations ([8]- 
ITT]) retaining only the terms corresponding to the central and first sidebands and obtain 
analytical results for central and first sidebands, to and t±i : 



tn 



cos 9o cos 00 cos ct] + i sin [A;^ a] ( 1 + sin 9o sin i 



2 Jo{a) cos (pn 

where n = ±1, t^o and tgn are transmission amplitudes for the static barrier at energy E 
and E + nhuj and 

.± _ A ± A ± 



A„ = cos[A;2 0. ± 6„]/ cos^„, 
Qn = isin[k2a]/ cos6'„. 

In the high barrier limit, |^| ^ with the result 0o,6n 0, we obtain expressions for 
transmission probabilities for the central and the sidebands. For the central band 

cos^ (f)0 

-^0 ~ ar, n 1 ■ 2 1 ~ -^^0 (13) 

1 — cos^[/e2aJ sm 0o 

where Tgo denotes the transmission probability at incident energy E and incident angle 0o 



in the case of the static barrier. This is the result obtained as Eq.(4) in HI]. For sidebands, 
we obtain: 

= cos[A;^a] ^ = -2 sin[(A;^ + A;°)a/2] sin[(A;^ - k°)a/2] 
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Qn = i sm[k^a] =^ A„ = 2i cos[(A;^ + k'^)a/2] sm[{k^ - k'^)a/2] 

tn = 2z4^^ sm[{k^ - k',)a/2]{cos[{k^ + ^2°)a/2] cos[(0„ - (/>o)/2] 

Jo{a) cos0„ 

+isin[(A;^ + A;°)a/2] cos[(0„ + (/,o)/2])e*('=?"+(^°-^")/2) 
The transmission probability for the sidebands is given by 

_ COs(0n) I |2 
COS(0o) 

V Jo(Q;j/ cos 0n cos 00 
+ sm'[{k^ + fc°)a/2] cos2[(0„ + 0o)/2]) 

where /k*; < £'cos0o otherwise TLi = 0. Tgn — \tsn\ is the transmission probabihty of 
electrons at energy E + nfuj and incident angle 0„ for the static barrier. We can also write 
the above expression as 

Tn = T^oTsn f24^y ^^^3^|^Ml^(sin0osin^ 
\ Jo[a) J cos 0„ cos 00 

(14) 

At normal incidence, 

7^±i=(2^)'sin^p2°-^r^)a/2] 
and ii fuj < \E — V\ we can write 

k^ - ki^ ^\E-V\ /hvF -\E-V±nu;\ /hvF = ±u;/vf 

with the result 

2 



2 

2 



r.:=|2%i^^l sin 

Jo{a) 



2'-fm si.n-/2] 

Jo(q;) J 



where t = a/vp is the time taken by a normally incident electron to cross the barrier without 
multiple refl,ections inside it. From the above expression, we note that Ti = T_i. For small 
a, J±i{a) ~ ±a/2; Jo{a) ~ 1 and sin[a;T/2] fa ujt/2 when cut is small, corresponding to low 
frequency limit where frequency is smaller than the reciprocal of the traversal time. Using 
these results we obtain 
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The above result can be compared with Eq.(8) in[4], where the transmission probabihty 
through a time-modulated barrier for the first sidebands is determined. The factor T, the 
transmission probability of the central band, is not unity and hence it appears there whereas 
TgQ = Tgn = 1, for normal incidence, in our case. 



A. Results and Discussions 



The results for the transmission of Dirac electrons in graphene are now presented. The 
following parameters were used: The Fermi wavelength of the incident electron is taken to 
be A = SOnm, the barrier oscillation frequency = 5 x the barrier width a = lOOnm 

and the barrier height V = 200meV. The dependence of transmission probabilities on a = 
Vi /hw for normally incident electrons and for those arriving at incident angle 30 degrees is 
shown in Figure(la,b), respectively. For normal incidence, the angular dependence of the 
transmission probability for the nth sideband is independent of the sign of n: T^n = T^n 
for fc]"'"' real. But this does not hold for incidence other than normal. We also find that 
the quantity a is very significant in determining the relative transmission probabilities of 
various sidebands as shown in the figure. This implies that by adjusting the value of a we 
can increase transmission through a particular sideband. It is seen that the central band 
dominates the transmission at all incident angles for small values of a whereas contributions 
from higher and lower sidebands increases as a becomes larger. This is plausible because for 
lower values of a the oscillating barrier can be treated as a static one since we are keeping 
uj fixed and changing Vi with the result that a is proportional to Vi in these figures. 
Moreover, the total transmission probability through the central as well as the sidebands is 
unity. Hence, perfect transmission for the oscillating barrier at normal incidence which was 
earlier observed for the static barrier This is due to the chiral nature of the particles 



which results in perfect transmission (Klein tunneling). 

In Figure (2a) we present the angular dependence of the transmission probability for the 
central-band Tq for various values of a = Vi /hw. The transmission probability for the static 
barrier is also shown in the figure as it corresponds to a = 0. The transmission probability T 
for the static barrier was previously obtained in|ll|. We find resonant transmission through 
the oscillating barrier but unlike the static barrier we do not find perfect transmission for 
any incident angle. Realize that for the static barrier there is perfect transmission for certain 
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values of the incident angle. This is to be expected as the probabilities are now spread over 
the central and sidebands. In addition, the maximum transmission through the oscillating 
barrier depends on the value of a. 

Figure (2b) shows the transmission probabilities for the central band along with the first 
few sidebands as a function of the incident angle for a = //lo; = 5. In this figure, we 
show how the incident particle flux is distributed in the sidebands (through the respective 
transmission probabilities) as the incident angle is varied. Note that the propagation angle 
for nth sideband is 0„ which is not the same as the incident angle 0o • For this particular 
value of a, transmission probability in the central band is very small for normal and close to 
normal incidence. For higher sidebands, more and more peaks in transmission probabilities 
occur. In the static case, the peaks in the transmission probability of the central band 
(there are no sidebands there) correspond to perfect transmission and the incident angles at 
which these occur can be obtained from the resonance condition, ki^ = ^ ( p is an integer), 
through Eq.( fT3l) and For the time- dependent situation being investigated here, it is not 
easy to determine the positions of the peaks as the analytic expression is more complicated. 
Nevertheless, we can understand how and where they occur by examining Eq ( |T4|) . albeit 
for small a where analytical results can be obtained but essential physics is the same. 
We observe, the transmission probability T„ given by Eq. (fl^ depends most strongly on the 
prefactor Tg^Tgn for the parameters considered here. The peaks correspond to the peak values 
of TsoTsn- Furthermore, the same behavior is seen for the static case as the transmission at 
higher incident energy there corresponds to transmission in the sidebands here. At these 
higher energies, the x-component of momentum in the barrier region satisfies the resonance 
condition greater number of times as the incident angle is varied, thus larger number of 
peaks. 

We note that the absence of any potential gradient along the ?/-direction results in the 
conservation of the y-component of momentum. Therefore, change in energy that an electron 
experiences due to exchange of modulation quanta with the oscillating barrier brings about 
corresponding changes only in the x-component of the electron's momentum. For non- 
zero fcy, energy exchanges can makes x-component of momentum imaginary inside or/and 
outside the barrier region that corresponds to unavailability of any energy state in the 
relevant region(s). If energy E + Ifkjj in the Ith sideband is such that \E + lhjjj\ < hvpky, 
there are no propagating states available outside the barrier since k{ becomes imaginary. At 
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the same energy when particles have states available inside the barrier it can be localized if 
it is transferred to these states after losing energy through interaction with the oscillating 
barrier. In this situation, the particles are confined across the barrier while they are free 
to propagate along the barrier till one or more quantum of energy is absorbed, allowing 
transition to a higher sideband with states aligned in energy outside the barrier leading to 
eventual escape from the barrier region. For a graphene quantum well, confined electron 
states which arise due to the suppression of electron-hole conversion at the barrier have been 
discussed in |l5|. 



For electron energy such that \E + Ihw — V\ < hvpky, there are no propagating states 
available inside the barrier since k2 becomes imaginary. Furthermore, the energy at which 
electronic states outside the barrier match the hole states inside it, electronic transmission 
is governed by Klein tunneling while unavailability of hole states inside the barrier results 
in ordinary tunneling. 

In Figure (3a) we present the transmission probability as a function of barrier width a 
for normal incidence. For the static barrier there is perfect transmission as can be seen in 
Figure (3a) where T represents the transmission probability for the static barrier whereas 
the transmission probability for the central band in the oscillating barrier decreases for 
smaller values of the barrier width and shows oscillatory but damped behavior for larger 
barrier width. The transmission probability for the other sidebands increases initially from 
zero but then oscillates with damped amplitude. We also observe that the contribution 
in transmission of the higher sidebands rises as the barrier width increases, this occurs 
due to larger time available to the electron for interacting with the oscillating potential as 
it traverses the barrier. In addition, we find that for normal incidence in the oscillating 
barrier: T+„ = T_„ for real. Nevertheless, the total transmission probability through 
the central as well as the sidebands is unity. These results imply that perfect transmission 
at normal incidence is independent of the barrier width, yet another manifestation of Klein 
tunneling. 

In Figure(3b), the transmission probability as a function of barrier width a when the 
incident angle is 30 degrees is shown. The transmission probability represented by T for 
the static barrier now oscillates as a function of the barrier width whereas transmission 
probabilities for the central and sidebands in the oscillating barrier show behavior close to 
that obtained for normal incidence. 
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A comparison between analytical result obtained in Eq (fT^ and numerical results is pre- 
sented in Figure (4) for a = 0.5. Transmission probabilities T±i of first sidebands are plotted 
against incident angle 0o- Inset shows plot of T_i versus 0o- K shows that transmission 
probabilities determined numerically exhibit the same behavior as obtained in the analyti- 
cal result. 

To summarize, we have considered the tunneling of chiral massless electrons corresponding 
to monolayer graphene through a barrier that is oscillating harmonically in time. We have 
determined how the transmission probability for the central and sidebands depends on the 
incident angle of the particles, the width of the barrier, the height and frequency with 
which it oscillates. Due to the chiral nature of the particles in graphene, tunneling is highly 
anisotropic with peculiar behavior at normal and close to normal incidence(Klein tunneling). 
We find, for normal incidence, perfect transmission in monolayer graphene. Klein tunneling 
that was observed for the static barrier is found to persist for the oscillating barrier. 
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